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ABSTRACT 

The ionization rate of interstellar material by cosmic rays has been a major source of controversy, 
with different estimates varying by three orders of magnitude. Observational constraints of this 
rate have all depended on analyzing the chemistry of various molecules that are produced following 
cosmic-ray ionization, and in many cases these analyses contain significant uncertainties. Even in the 
simplest case (H^~) the derived ionization rate depends on an (uncertain) estimate of the absorption 
path length. In this paper we examine the feasibility of inferring the cosmic-ray ionization rate using 
the 10830 A absorption line of metastable helium. Observations through the diffuse clouds toward 
HD 183143 are presented, but yield only an upper limit on the metastable helium column density. A 
thorough investigation of Hc + chemistry reveals that only a small fraction of He + will recombine into 
the triplet state and populate the metastable level. In addition, excitation to the triplet manifold of 
helium by secondary electrons must be accounted for as it is the dominant mechanism which produces 
He* in some environments. Incorporating these various formation and destruction pathways, we derive 
new equations for the steady state abundance of metastable helium. Using these equations in concert 
with our observations, we find £hc < L.2 x 10~ 15 s _1 , an upper limit about 5 times larger than the 
ionization rate previously inferred for this sight line using H3 . While observations of interstellar He* 
are extremely difficult at present, and the background chemistry is not nearly as simple as previously 
thought, potential future observations of metastable helium would provide an independent check 
on the cosmic-ray ionization rate derived from H^ in diffuse molecular clouds, and, perhaps more 
importantly, allow the first direct measurements of the ionization rate in diffuse atomic clouds. 
Subject headings: astrochemistry — atomic processes — cosmic rays 



1. INTRODUCTION 

1.1. Motivation 

Over the past several decades, the assumed value of 
the cosmic-ray ionization rate of interstellar hydrogen 
has fluctuated up and down. Various theories and 
models have predicted ionization rates from 10~ 18 s _1 
to 10~ 15 s -1 in the diffuse interstellar medium (e.g . 
Spitzer fc Tomaskol 119681: Ivan Dishoeck fc Blackl 119861 : 



Webber! 1 19981 iLe Petit et alj |2004| )T On the other 
hand, observations of molecules such as HD and OH 
typically resulted in estimates of t he ionization rate that 
were on the order of 10~ 17 s" 1 (lO'Donnell fc Watson! 



1971 iBlack fc Dakarnd H977t Black et al l 11971 



Hartquist et alJ Il978al |bl: iFederman et al.l 119961 ). How- 



ever, these estimates depend on gas phase abundances 
of O, OH, D, and HD, values which are often difficult 
to measure precisely. More recently, observations of 
Hg" have again revised the cosmic-ray ionization rate 
upward to a few tim es 10" 16 s" 1 (jMcCall et al.1 120031 : 
Ilndriolo et al.l l2007f ). Deriving the ionization rate 
from H^ requires only one uncertain parameter, the 
absorption path length. While the higher ionization 
rates derived fro m H;j" are becoming generally accepted 
(Dal garno I l2006h . it is desirable to search for new ob- 
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servables which can offer independent and less uncertain 
estimates of the ionization rate. 

In this report, we investigate the possibility of obser- 
vationally determining the total ionization rate of helium 
atoms by cosmic rays in diffuse clouds. The basic premise 
is that in a sufficiently reddened cloud, the column den- 
sity of neutral helium atoms excited to the metastable 
ls2s 3 Si level may be high enough to be measured by 
means of interstellar absorption lines arising from that 
level. The high cosmic abundance of helium and the 
long radiative lifetime of the metastable level, A^ 1 = 2.5 
hr, may compensate for the difficulty of populating this 
highly excited level, which lies 19.8 eV above the ground 
level. Previously, it has been assumed that this level 
should be populated almost entirely by cosmic-ray ion- 
ization of helium atoms, followed by radiative recombi- 
nation of the ions with electrons. Figure Q] schematically 
shows the processes conventionally used in describing the 
(de)population of the metastable state. 

1.2. Background 

A simple reaction network - consisting of (i) cosmic- 
ray ionization of He atoms in the Is 2 ground level, 
(ii) radiative recombination of He + ions with electrons to 
the metastable ls2s 3 Si level, and (iii) forbidden spon- 
taneou s emiss ion to the ground level - was first used by 
iScherbl (p6Sh and iRees et al.l ljl968h in proposing the 
observability of interstellar metastable helium. However, 
both of these studies considered 2-photon emission from 
the metastable state, then thought to be the dominant 
method of depopulation. This 2-pho ton decay proceeds 
on a timescale of A^ 1 > 116 days (Drake fc Dalgarnol 
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Fig. 1. — Simplified energy level diagram of atomic helium, 
showing only S, P, and D terms up through a principal quantum 
number of 10. The key processes thought to control the abundance 
of the metastable 2 3 Si state, along with the absorption line from 
this state at 10830 A, are indicated. It is assumed that all electron 
recombinations into either the singlet or triplet manifold quickly 
lead to the ground state of that manifold by allowed spontaneous 

¥§6cF)° n ' Iii the following year, the radiative lifetime 
associated with 1-photon dec ay was computed to be 
A^ 1 w 7 hours (|Grieml fl969f ). nearly 400 times faster 
than the 2-photon decay. This value was later re- 
fined to A' 1 = 2.5 hours (jWoodworth fc Moos! Il975t 
lHata fe GrantJ[l98U ). still much faster and as a result 
the an alyses performed bv lScherbl (|1968fl and lRees et al.1 
(1968) had greatly overestimated the popula tion in the 
meta stable state (this possibility was noted bv lRees et al.l 
1968). Consequently, the thought of observing interstel- 
lar metastable helium was abandoned. 

Because of the high ionization rate inferred from Hjj" , 
we decided to revisit these calculations considering up- 
to-date rate coefficients and improved telescope/detector 
capabilities. Assuming the same chemical scheme as in 
the past, we can derive the steady state equations for the 
ground, ionized, and metastable states: 



ain e m + n m A = n g ( He , 
rigCiie = («i + a 3 )n e ni, 
a 3 n e m = n m A, 



(1) 
(2) 
(3) 



Here, n g , rij, and n m denote the populations of the 
ground, ionized, and metastable levels, respectively; ^He, 
the total ionization rate of He atoms due to cosmic rays, 
including the effects of secondary electrons; ai and 0:3, 
the total, direct recombination rates to all singlet levels 
and to all triplet levels, respectively; n e , the electron den- 
sity in the gas; and A, the Einstein coefficient for spon- 
taneous emission from the metastable level. The units 
of each term in equations |T]) through ([3]) are cm~ 3 s -1 . 
All recombinations to triplet levels above the metastable 
level are assumed to produce subsequent radiative cas- 
cades to the metastable level that are effectively instan- 
taneous, owing to the long lifetime of the latter level. 



Similarly, all recombinations to singlet levels are assumed 
to cascade promptly to the ground level. Equations 
and ([3]) can be solved for the ratios ni/n g and n m /ni, 
and thus for n m /n g as well. These ratios can be substi- 
tuted into the definition of the fractional population of 
the metastable level, 



(4) 
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in order to obtain the desired relation between the frac- 
tional metastable population f m and the ionization rate 

CHo, 



J m 
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(a 3 n e ) (>Chc) ' 



(5) 



The triplet branching fraction for recombinations at 
70 K is b — 0:3/(0:1 + a 3 ) — 0.62, since ct\ — 4.0 x 
10~ 12 cm 3 s" 1 and a 3 = 6.6 x 10~ 12 cm 3 s" 1 (R. Porter 
2009, private communication). Radiative decay to the 
ground level is by far the fastest of the three processes 
mentioned above, with A = 1.1 x 10 -4 s _1 . In con- 
trast, n e = 0.02 cm -3 , 03n e = 1.3 x 10~ 13 s _1 , and 
Chc = 3 x 10~ 16 s _1 are representative values in diffuse 
clouds. Given these values, 1 <C A/(a 3 n e ) <C A/(b(n c ), 
and equation {5} can be approximated by 
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Owing to the very large differences among the rates, this 
approximation to equation ([5]) is nearly exact. This holds 
true as long as £hc <C n e (ai+a 3 ) ~ 10~ 13 s _1 , such that 
ionization of helium by cosmic-rays is the rate-limiting 
step in the path to the metastable state. In this limit 
f m effectively depends on £hc alone - apart from the 
well-determined atomic constants b and A - thus sug- 
gesting metastable helium as a fairly robust indicator of 
the cosmic-ray ionization rate. 

2. OBSERVATIONS 

2.1. Predictions 

The fundamental question remaining then is whether 
the interstellar lines of He* arising from a suitable dif- 
fuse cloud are likely to be detectable. The strengths of 
these lines are fixed by the cloud's column density of 
metastable atoms, N m , which can be calculated from 



N m = f m N(He) = f m A(He)Nn, 



(7) 



where iV(He) is the total column density of helium atoms 
in all states, iVn is the total column density of hy- 
drogen nuclei [JV H = A(H) + 2jV(H 2 )], and A (He) = 
N(He)/N-a — 0.097 is the re lative abundance of helium 
with respect to hydrogen (|Anders fc Grevessi 119891) . 
The fraction of interstellar helium sequestered in the 
grains has also been assumed negligible. If a direct mea- 
surement of Nn is not available, an alternative is to use 
Ah = (3E(B — V), where E{B — V) is the observed color 
excess, and /3 = N H /E{B-V) = 5. 8 x 10 21 cm" 2 mag" 1 
is the interstellar gas-to-dust ratio (|Bohlin et al.lll978T ). 

To estimate the expected line strengths, we assume 
E(B - V) = 1.0 mag and Ch = 3 x 10" 16 s" 1 in a suit- 
able, individual interstellar cloud. The former value leads 
to Nk = 5.8 x 10 21 cm" 2 and N(He) = 5.6 x 10 20 cm' 2 . 
A substitution of the assumed value of Chg into equation 
© gives f m - 1.7 x 10- 12 . Then, N m = f m N{Re) = 
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TABLE 1 
The 1s2s 3 S - 1s2p 3 P Multiplet of 
He i 











(A) 


J(lowcr) 


J (upper) 


/ 


10829.0911 


1 





0.060 


10830.2501 


1 


1 


0.180 


10830.3398 


1 


2 


0.300 



Note. - Wavelengths and oscillator 
strengths are from the NIST Atomic Spec- 
tra Database (Ralchenko ct al. 2008) 



9.5 x 10 8 cm~ 2 . The best choice among the available 
He I* lines is the ls2s 3 S - ls2p 3 P multiplet located near 
10830 A. Data for the transitions associated with this 
multiplet are shown in Table 1 , where column 4 gives the 
oscillator strengths. These lines are stronger than other 
transitions arising from the metastable level (such as the 
multiplet near 3889 A), and the near-infrared wavelength 
is advantageous in observations of heavily reddened stars 
with large total column densities of helium. 

Assuming N m = 9.5 x 10 s cm~ 2 , the equivalent width 
of an unresolved blend of the two strongest lines of the 
multiplet, which are separated by only 2.5 km s _1 , would 
be W\ = 0.47 mA. If a spectrometer with a resolving 
power of 70,000 were used, the line would have a central 
depth of ^0.30%, thus demanding a signal-to-noise ratio 
(S/N) of ~1000 on the continuum for a 3cr detection. 
Modern optical echelle spectrographs can easily reach 
S/N exceeding 2000 (e.g. JAdamkovics. Blake, fc McCaiH 
2003), but reaching such a high S/N in the near-infrared 
is a significant challenge. 

2.2. Target Selection 

In choosing a target, we searched for sight lines that 
had a combination of several desirable characteristics: 
high color excess; high cosmic-ray ionization rate as in- 
ferred from H^"; relatively bright J-band magnitude; few 
interstellar velocity components; well-behaved stellar ab- 
sorption features. Using these criteria, we arrived at 
HD 183143 as our most favorable target, with J = 4.18, 
V = 6.86, E(B - V) = 1.27, and a spectral type of 
B7Iae. The star's photospheric He I absorption lines 
at 5875, 6678, and 7065 A are relatively broad, with 
FWHM > 66 km s _1 , and they are also free from an y 
interfering emission components (|Thorburn et al.1 l2003). 
This suggests that the photospheric 10830 A line may 
desirably provide a smooth background with a shallow 
slope, against which one could search for the much nar- 
rower, weak interstellar lines. In addition, many inter- 
stellar absorption lines along this line of si ght have been 
previously studied, including CH and CH+ (jGredel et al.l 
[l99"3h . and CN, H+, and 12 CO (jMcOall et all [2002]) . 
These observations reveal that there are two distinct in- 
terstellar cloud groups at different velocities. Although 
this means that not all of the interstellar helium is at one 
velocity, it does provide a very useful method for poten- 
tially confirming a detection. Additionally, H;j~ observa- 
tions have been used to determine the ion ization rate of 
mole cular hydrogen, £2, in this sight line (In driolo et al.l 
12001 . 

Using values specific to the HD 183143 sight line 



(E(B -V) = 1.27; Ch c = 3.5 x 10~ 16 s" 1 ) we can again 
perform the calculations in §2.1 to determine the ex- 
pected line strength. The resulting equivalent width is 
W\ = 0.70 mA. Because the sight line has 2 velocity 
components though, we assume equal amounts of ma- 
terial in each cloud, and thus expect 2 absorption lines 
with W\ = 0.35 mA. These would require S/N~1300 
for a 3cr detection given the same instrument capabilities 
assumed above. While obtaining a S/N this high is diffi- 
cult in the near-infrared, some of the most advanced tele- 
scope/detector combinations are capable of approaching 
such results, so we proceeded with observations. 

2.3. Execution 

Data were taken near the He I* line at 10830 A us- 
ing the Phoenix spectrometer (jHinkle et a l. 2002) on the 
Gemini South Telescope. The spectrometer was used 
with its echelle grating and 0.17" slit to produce a re- 
solving power of ^70,000, and the J9232 filter to select 
the correct order. Observations of both the target (HD 
183143) and standard (a Aql) stars were taken on May 
25, 2008 and June 28, 2008. Total integration times were 
33 and 30 minutes for the target and 1.9 and 1.4 min- 
utes for the standard on each night, respectively. During 
each set of observations, the star was nodded along the 
slit in an ABBA pattern to allow for the later subtrac- 
tion of neighboring images, and thus the removal of the 
atmospheric background and detector bias levels. 

3. DATA REDUCTION 

A significant portion of the data reduction - dark cur- 
rent subtraction, subtraction of neighboring images, re- 
moval of bad pixels, flat fielding, combination of ex- 
posures with the spectral image in the same nod posi- 
tion, fitting of the spectral response, and spectral ex- 
traction - was performed using NO AO's IRAF package 5 . 
Du ring thi s proce ss, we combined the methods outlined 
by iKulesal (|2002l ) with those given by NO AO's online 
Phoenix documentation 6 in order to obtain the best pos- 
sible S/N. Once the one-dimensional spectra were ex- 
tracted, they were imported to IGOR Pro 7 , w here we 
have macros set up to complete the reduction (|McCalll 
l200lh . 

Because of the annual shift in the relative positions 
of (inter)stellar and atmospheric features with time, the 
data from different nights were first analyzed separately. 
In all cases, however, the expected locations of the inter- 
stellar He I* lines lie within the broad stellar absorption 
line. Because the S/N of the standard star was actually 
lower than that of the target in the June data, we de- 
cided to forego the process of ratioing by the standard 
star and we instead directly fit the stellar absorption fea- 
ture. Both the A and B beams for each night were wave- 
length calibrated using atmospheric lines and then added 
together. The broad stellar absorption feature was then 
fit using the summation of three gaussian functions, all of 
which were constrained to have FWHM at least 3 times 
that of the 10 km s _1 measured for interstellar absorp- 
tion features along the line of sight. The spectra from 
each night were then divided by their respective fits and 

5 See http://iraf.noao.edu/! 

6 See http://www.noao.edu/usgp/phoenix/phoenix.html 

7 See http://www.wavemetrics.com/ 
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Fig. 2. — Spectra of HD 183143 in velocity space. The top spec- 
trum, observed with Phoenix at Gemini South, has been adjusted 
for the centroid of the unresolved blend of the 1-1 and 1-2 mem- 
bers of the He I* multiplet, and the broad p hotospheric line has 
been divided out. The bottom spectrum (from.McCall ct al1 l2002r) 
shows the R(l, l) 1 transition of H^" for reference. Arrows indicate 
the two interstellar velocity components which have been observed 
in various molecules (CH, CH+, CN, 12 CO, H^~). 

shifted to be in the local standard of rest (LSR) frame. 
Finally, the fully reduced spectra from both nights were 
added together and converted to velocity space to pro- 
duce the top spectrum shown in Figured 

4. RESULTS 

There is no indication of interstellar He* absorption at 
either of the expected velocities. While we did obtain a 
relatively high S/N (~ 700) for high-resolution infrared 
spectroscopy, we were unable to achieve the desired S/N 
~ 1300. The non-detection of the He I* lines enabled 
us to calculate an upper limit to the column density of 
metastable helium along this line of sight. First, the 
upper limit to the equivalent width, W\, was computed 



(8) 



where a 
trum, A 



pix 



0.00145 is the standard deviation in the spec- 
i,,x = 0.05 A is the wavelength per pixel, and 
13 is the number of pixels expected in a single 
absorption component given a 10 km s _1 FWHM (this is 
the average FWHM of H+ lines reported bv lMcCall etldl 
(2002)). These quantities result in an upper limit to the 
equivalent width of W\ < 0.26 mA (W\ < 0.78 mA at 
the 3er level). 

Next, the standard relation between equivalent width 
and column density was used: 



N = 



W x m e c 2 
A 2 7re 2 / 



(9) 



where m e is the electron mass, c is the speed of light, A is 
the wavelength of the transition, e is the electron charge, 
and / = 0.48 is the sum of the oscillator strengths of the 2 
strongest, blended lines. Inserting the known parameters 
results in a 3cr upper limit to the column density AT m < 
1.6 x 10 9 cm -2 in a single velocity component. However, 
there are 2 cloud components along this sight line, so the 
total line of sight limit is N m < 3.2 x 10 9 cm~ 2 . 

5. ANALYSIS 
5.1. Reaction Network Revisited 



Hc+ + H -> Hc + H+, 
Hc+ + H 2 -> Hc + H+, 



In planning observations and calculating predicted 
line strengths, we relied on the simple chemical scheme 
that only considers the destruction of He + via elec- 
tron recombination. However, during the course of this 
study we ident i fied (f rom advanced chemical models, e.g. 
iWoodall et al.1 (|2007t )1 several competing reactions that 
could be important in destroying He + in diffuse clouds: 

(10) 
(11) 

H 2 ^Hc + H + H+, (12) 
CO->He+0 + C+. (13) 

Rate coefficients for these reactions, as well as the elec- 
tron recombination reactions, can be determined for a 
specific temperature, T (in Kelvin), using the fitting pa- 
rameters (a, /3, and 7) in Table in conjunction with 
the expression 

* ( T Y 



He+ 
He+ 



cm 3 s" 



(14) 



Unlike the case of electron recombination, these re- 
actions should not lead to metastable helium. Simple 
energetics arguments demonstrate why this is the case. 
The energy difference between the ionization potential 
of helium (24.6 eV) and the excitation energy of the 
metastable state (19.8 eV) is only 4.8 eV. In order to 
dissociate and/or ionize the reaction partners of He + , 
reactions (fTUl - [T3"|) require energies of 13.6, 15.4, 18.1, 
and 22.5 eV, respectively (assuming all reactants are in 
the ground electronic state). At diffuse cloud tempera- 
tures (~70 K) thermal energy is much less than 1 eV, 
and so the kinetic energy of the reactants will have no 
effect. Following these arguments, the neutral helium 
product can only be in the ground state as more than 
4.8 eV of the helium ionization potential must be used in 
each reaction. As a result, accounting for these reactions 
greatly decreases the fraction of helium ions which pass 
through the triplet manifold, and thus the population of 
the metastable state. 

However, we also made the assumption that the 
metastable state is only populated via cosmic-ray ion- 
ization of He, followed by electron recombination. Given 
that a smaller branching fraction limits this pathway, 
electron impact excitation into the triplet manifold will 
be a competing formation mechanism. Cross sections for 
ionization and excitation of heli um by electrons in the 
10 — 1000 eV range are shown in iDalgarno et ail (|1999f ) 
figures 2a & 2b. To compute the rate of ionization and 
excitation, one must perform an integral in energy space 
over the product of each cross section with the differen- 
tial energy spectrum of electrons in the ISM. This full 
calculation is hindered by the fact that the spectrum 
of secondary electrons (those produced during ionization 
events) is unknown, and cannot be derived from the dif- 
ferential energy spectrum of cosmic-ray protons which is 
also unknown below ~ 1 GeV. The complexity associ- 
ated with deriving the spectrum of secondary electrons 
is beyond the scope of this paper, thus we make some 
simplifications in estimating the importance of electron 
impact excitation into the triplet manifold of helium. 

Assuming that all secondary electrons have the same 
energy, the ratio between the rate of excitation into all 
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triplet states and the rate of ionization can be deter- 
mined by taking the ratio of the respective cross sections 
at a given ene rgy. We take this rati o at 30 eV (the mean 
value given bv lCravens fc Dakarnol (|1978f) 1. and find the 
rate of excitation into all triplet states to be 2 times faster 
than the rate of ionization by secondary electrons. To de- 
termine the overall importance of electron impact excita- 
tion then, we need to find a relationship between the total 
ionization rate of helium and the ionization rate due to 
secondaries. Using relations betwe en the primary ioniza- 
tion rates of hydr ogen and helium (|Habing fc Goldsmith! 
119711 iLisztl 12003ft and between the primary ionization 
rate of hydrogen and the t o tal ionization rate of helium 
(iGlasseold fc Laneeri Il97i iTielend [2005), we estimate 
that ionization by secondary electrons accounts for about 
1/6 of the total ionization rate of helium. This, in turn, 
leads to the approximation that the rate for electron im- 
pact excitation into the triplet manifold — and thus the 
metastable state (which we denote fee*) — should be 
roughly 1/3 that of the total ionization rate of helium 
(i.e. fe e * ~ Cho/3; we use this relation for the remainder 
of this paper). 

Mathematically, these additional formation and de- 
struction reactions can easily be included by altering the 
steady state equations in §1.2, resulting in 2 changes to 
our analysis. First, due to the additional destruction 
pathways of He + , the branching fraction, 6, must be re- 
defined as 

, _ O^e 

~ n(H)^+n(H2)(^+^l2i)+n(CO)^l3]+n e (a 1 +a3)■ 

(15) 

In many cases, however, absolute abundances are not 
known and it is thus convenient to recast equation (|15jl 
in terms of fractional abundances as 

b _ azx e 

( i - /h 2 ) tan+ /h 2 (*rrn+ ¥rj) / 2 + ^co *rrj+ x e («i + « 3 ) ' 

(16) 

where Xj = nj/nn, riH = u(H) + 2n(H2), and the molec- 
ular hydrogen fraction /h 2 = 2n(fi2)/riH- Second, equa- 
tion (J6j) must be recast to include electron impact exci- 
tation into the metastable state, and becomes 



fm Kttc + fee* 

While the analysis now includes many more parame- 
ters, we can still calculate the fractional abundance of 
metastable helium, and thus the expected line strength, 
toward HD 183143. We assume that fractional abun- 
dances are constant throughout the cloud, allowing us to 
substitute column densities for number densities when 
available (i.e. Xj = Nj/Nn). Using the color excess 
as in §2.1 gives ife = 7.4 x 10 21 cm~ 2 . This is used 
in conjunction with spectroscopic observations of CO 
which indicate N(CO) « 10 15 cm" 2 (jMcCall et alj|2053) 
to compute xco- The assumption that there are equal 
amounts of atomic and molecular hydrogen is quanti- 
fied by /h 2 = 2/3. Finally, observations of C + in dif- 
fuse clouds have shown that x e ~ 1.4 x 10 ~ 4 , assum- 
ing that nea rly all electrons com e from this singly ion- 
ized carbon ( Cardelli et al.l lT996) . Combining these data 
and assumptions with the rate coefficients in Table [H 
the new branching fraction is b = 0.08, about one- 
eighth of the value considering electrons alone. Sub- 



stituting this branching fraction and the relevant pa- 
rameters from §2.1-2.2 into equation (jTTJ) results in val- 
ues of f m = 1.3 x 10~ 12 , N m = 9.3 x 10 s cm" 2 , and 
W\ = 0.46 ml. A gain splitting the material into 2 
equal cloud components decreases the equivalent widths 
to W x = 0.23 ml, which would require a S/N~ 2000 for 
a 3<7 detection. 

5.2. Cosmic-Ray Ionization Rate of Helium 

Re-arranging equation (jTTJ) , we can turn this prob- 
lem around and compute an upper limit to the cosmic- 
ray ionization rate of helium using our observations. 
Given the upper limit to the metastable column density, 
N m < 3.2 x 10 9 cm~ 2 , and the estimated total helium 
column, AT(He) = (3A(He)E(B - V) = 7.1 x 10 20 cm" 2 , 
the 3cr upper limit to the fractional metastable pop- 
ulation is f m < 4.5 x 10~ 12 . Using this in concert 
with the branching fraction above, b — 0.08, results in 
fe < 1.2 x 10 -15 s _1 . This upper limit is about 5 
times larger than the ionization rate inferred from H^ 
observations (assuming the relation between the ioniza- 
tion rate of heliu m and molecular hydrogen is given by 
2.3C He = 1.5Ca (|Glasseold fc Langerl fl97l ). Because 
of electron impact excitation into the metastable state 
though, this determination of the ionization rate relies 
on a much more indirect analysis than was initially pro- 
posed. Comparing &fe c to fe c *, we can compute the 
relative importance of each formation mechanism via 

= TTr 6 Tx = (3^ + (18) 

Ofec + OHe* 

In doing so, we find that electron impact excitation ac- 
counts for 80% of the metastable population while ion- 
ization and electron recombination accounts for 20%. 

6. DISCUSSION 

While the reactions associated with metastable helium 
are more complex than previously presented, we still see 
it as a viable tracer of the cosmic-ray ionization rate. 
As such, we decided to investigate the prospects for He* 
detections in various interstellar environments, includ- 
ing diffuse molecular clouds (100 < jih S 500 cm~ 3 , 
/h 2 0.1), dense clouds (tih > 10 4 cm~ 3 , /h 2 ~ 1), 
and diffuse atomic clo uds («h ^ 100 cm" 3 , /h 2 < 0.1) 
(jSnow fc McCal1ll2006fh The following analyses will high- 
light the branching fraction in each environment, as well 
as the relative importance of electron impact excitation 
vs. ionization + recombination using equation (I18|) . 

6.1. Diffuse Molecular Clouds 

Given that the analysis in §2.1 did not account for 
the processes examined in §5.1, we felt it prudent to re- 
visit the calculations for diffuse molecular clouds. We 
use the same values as before (E(B — V) = 1; fee = 
3 x 10~ 16 s _1 ), but now also assume /h 2 = 2/3, x e = 
1.4 x 10~ 4 , and xco = 10~ 7 . The general results for 
this environment (b = 0.08; f m = 1.1 x 10~ 12 ; N m — 
6.3 x 10 s cm -2 ; W\ — 0.31 mA) are similar to those 
for the specific diffuse molecular sight line HD 183143, 
with the differences due to the lower color excess. As- 
suming that all of the material has the same velocity, 
metastable helium absorption should be observable in 
diffuse molecular clouds at a 3<7 level with S/N~ 1500. 
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Given the small branching fraction above, P(<5He*) = 0.8 
and we conclude that metastable helium is predomi- 
nantly formed via electron impact excitation in diffuse 
molecular clouds. 

6.2. Dense Clouds 

Dense clouds, while providing a larger total helium 
column, have several characteristics detrimental to the 
formation of metastable helium. The cosmic-ray ioniza- 
tion rate tends to be about 1 order of magnitude lowe r 
in dense clouds than in diffuse clouds (Dalgarno 2006). 
Also, the fractional abundance of electrons is much lower, 
x e « 4 x 10~ 8 , while the fraction al abundance of CO i s 
much higher, x co » 1.4 x 10~ 4 ()Woodall et al.ll2007D . 
Because l f[$\ is so much larger than any of the other rate 
coefficients, collisions with CO will dominate the destruc- 
tion of He + and equation (fTo| can be simplified to 

6«^^. (19) 
zcofcco 

Given the fractional abundances above and the rele- 
vant rate coefficients («3 and /ceo were computed for 
T = 40 K), the branching fraction is b ~ 10 -6 . As a 
result, P((5hc*) ~ 1, meaning that metastable helium is 
formed exclusively by electron impact excitation in dense 
clouds. Even with this formation mechanism though, the 
expected equivalent width (W\ — 0.13 mA) and neces- 
sary S/N for a 3cr detection (S/N~ 3700), coupled with 
the large attenuation of the background star's flux at 
1 /urn, make the detection of He* in dense clouds highly 
unlikely. 

6.3. Diffuse Atomic Clouds 

Diffuse atomic clouds, on the ot her hand, have negligi- 
ble concentrations of H 2 and CO (jSnow fc McCall 2006) 
and presumably share the high ionization rate of diffuse 
molecular clouds. In purely atomic conditions, electron 
recombination only has to compete with reaction pop 
and equation (|1€>[) can be approximated as 

b~i V 5 \ ' ( 2 °) 

TTQI+ x e {ai + a 3 ) 

The simplified result for atomic clouds is then b « 0.40, 
with a corresponding P((5h * ) = 0.45, meaning that 
ionization and electron impact excitation play roughly 
equal roles in forming metastable helium in such envi- 
ronments. Despite this branching fraction being closer 
to the ideal case of b = 0.62, the low amount of mate- 
rial along such a sight line (E(B — V) ~ 0.1) results in 
a predicted equivalent width of W\ rs 0.06 mA. How- 
ever, there are some diffuse atomic sight lines with more 
favorable conditio ns. One such candidate, a Sco, has 
E{B - V) = 0.40 (|Clavton fc Hanson! 119931 ) and thus a 
predicted equivalent width of W\ « 0.22 mA using equa- 
tions ([17]) & ([20|) . However, a Sco also has measured val- 
ues of7V(H) = 2.4 x lO 21 cm" 2 , iV(H 2 ) = 6.2 x 10 19 cm" 2 
([Savage et alJ fl977h . and AT(CO) = 6.5 x 10 12 cm" 2 
(I Allen et al.l l 1990'). which correspond to /h 2 = 0.049 and 
xco = 2-6 x 10~ 9 . Using these values and equation (fl6|) . 
we can test the accuracy of equation (|2H)) at small molec- 
ular fractions. The result is b = 0.31, or about a 30% 
error in the approximation. At fn 2 = 0.15, equation 
(|2"0|) overestimates b by a factor of 2, so this approxima- 
tion should only be applied for /h 2 < 0.1. Taking the 



branching fraction from the full calculation, we predict 
an equivalent width of W\ « 0.20 mA, and a correspond- 
ing S/N ~ 2400 necessary for a 3cr detection. If such a 
detection can be made, however, it will provide the excit- 
ing opportunity to probe the cosmic-ray ionization rate in 
an environment where H^" observations cannot be made 
due to the low molecular fraction. 

7. CONCLUSIONS 

We have analyzed the possibility of detecting ab- 
sorption lines due to interstellar metastable helium at 
10830 A. Observations toward the diffuse cloud sight line 
HD 183143 were taken, and a spectrum with S/N~700 
was obtained, but no interstellar He 1* lines were de- 
tected. In examining the chemistry associated with 
metastable helium, we have identified important forma- 
tion and destruction pathways, and have derived new 
equations for the steady state analysis of He*. While 
these reactions have been known for some time, this is the 
first instance where they have been applied to metastable 
helium chemistry. Using our observations and the newly 
derived equations, we inferred an upper limit for the 
cosmic-ray ionization rate of helium which, although con- 
sistent with other studies, is about 5 times larger than 
previously inferred values. 

To determine if future observations of interstellar 
He* are warranted, we predicted the S/N ratios nec- 
essary for 3cr detections in various environments. Dif- 
fuse molecular clouds are the most promising targets 
with S/N~1500 required, while favorable diffuse atomic 
clouds need S/N~2400. While such observations are 
extremely challenging at present, advancements in tele- 
scope and near-infrared detector technology may one day 
make metastable helium a widely applicable probe of the 
cosmic-ray ionization rate. In diffuse molecular clouds, 
He* will act as a cosmic-ray probe independent of H3", 
and together with H^ it will also enable determination 
of the absorption path length and average cloud density. 
He* observations will also be especially important for 
diffuse atomic clouds, where there are no other reliable 
tracers of the ionization rate. 

The authors wish to thank Tom Kerr and the UKIRT 
Service Observing Programme for obtaining preliminary 
spectra near 10830 A using CGS4, Ryan Porter and Gary 
Ferland for supplying unpublished values of the recom- 
bination coefficients and for suggesting the possibility of 
electron impact excitation into the metastable state, and 
the anonymous referee for helpful suggestions and com- 
ments. NI and BJM have been supported by NSF grant 
PHY 05-55486. 

Based on observations obtained at the Gemini Obser- 
vatory, which is operated by the Association of Univer- 
sities for Research in Astronomy, Inc., under a cooper- 
ative agreement with the NSF on behalf of the Gemini 
partnership: the National Science Foundation (United 
States), the Science and Technology Facilities Coun- 
cil (United Kingdom), the National Research Council 
(Canada), CONICYT (Chile), the Australian Research 
Council (Australia), Ministerio da Ciencia e Tecnologia 
(Brazil) and Ministerio de Ciencia, Tecnologia e Inno- 
vation Productiva (Argentina). The Gemini/Phoenix 
spectra were obtained through program GS-2008A-Q-14. 

The observations were obtained with the Phoenix in- 
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frared spectrograph, which was developed by the Na- tional Optical Astronomy Observatory. 
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TABLE 2 

Rate Coefficients for Reactions Involving Ionized Helium 









Coefficient at 70 K 


Reaction 


a 




7 (cm 3 s -1 ) References 



He+ - 


h H - 


> He 4 


- H+ 




1.2 X 10~ 15 


0.25 





>m= 


8.3 x 10- 


-16 


1 


He+ - 


h H 2 - 


-> He 


+ H+ 




7.2 x icr 15 










7.2 x 10" 


-15 


2 


Hc+ - 


h H 2 - 


-+ He 


+ H - 


-H+ 


3.7 x 10 -15 





35 




2.2 x 10" 


-14 


2 


Hc+ - 


h CO 


— He 


+ O 


+ C+ 


1.6 x 10" 9 










= 1.6 x 10 


-9 


3.1 


He+ - 


h e -> 


He(l 


1 S) + 


hv 


1.76 x icr 12 


-0.56 





a\ = 


4.0 x 10" 


12 


5 


Hc+ - 


h e -> 


He(2 


3 S) + 


hu 


2.84 x 10~ 12 


-0.59 





« 3 = 


6.6 x 10" 


12 


5 



References. — (11 IStancil et all (1998T) : f2l lBarlowl lfl984T) ; (3) ILaudenslaeer et al] Jl97^) ; 
(4) lAnicich et al.l 1119771 ); (5) R. Porter 2009, private communication 

Note. — Coefficie nts a t temperatures between about 10 and 300 K can b e de ri ved usin g 
a, (3, 7, and equation \1A\ . Rate coefficients and their references for reactions HOD , mil l, 11121 . 
and 1131 were found at http://www.udfa.net/ 



